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ABSTRACT 

We report the discovery of a z ~ 9 Lyman Break Galaxy (LBG) candidate, selected from the NICMOS 
Parallel Imaging Survey as a J-dropout with /[ \q - H lb o =1.7. Spitzer/IRAC photometry reveals that the galaxy 
has a blue H\(,q - 3.6 fj,m color, and a spectral break between 3.6 and 4.5 /im. We interpret this break as the 
Balmer break, and derive a best-fit photometric redshift of z ~ 9. We use Monte Carlo simulations to test the 
significance of this photometric redshift, and show a 96% probability of z > 7. We estimate a lower limit to the 
comoving number density of such galaxies at z ~ 9 of 4> > 3.8 x 10~ 6 Mpc" 3 . If the high redshift of this galaxy 
is confirmed, this will indicate that the luminous end of the rest-frame UV luminosity function has not evolved 
substantially from z ~ 9 to z ~ 3. Still, some small degeneracy remains between this z ~ 9 model and models 
at z ~ 2 - 3; deep optical imaging (reaching 7^ ~ 29) can rule out the lower-z models. 

Subject headings: galaxies: high-redshift - galaxies: evolution - galaxies: formation 



1. INTRODUCTION 

The search for the first galaxies remains one of the primary 
goals of extragalactic astronomy. Surveys of distant star form- 
ing galaxies (e.g. Lyman break galaxies; LBGs) have resulted 
in determinations of lumi nosity functions, stellar masses, and 
ages o f objects to z ~ 6 (Stanwav et alJ l2003h iBunker et alJ 
20041: lYan eTaTJ l2006t fYoshida et all 120061: iBou wens et all 
20061 lEvles et all 120071: IStark et al.ll2007alh The space den- 
sity of galaxies at z ~ 6 and their apparent maturity suggests 
that galaxy formation must have begun at z > 7. Although 
the requirement of deep near-infrared imaging makes z > 7 
difficult to observe, this is a critical epoch in the history of 
galaxy formation. Discovery of the first galaxy-wide bursts of 
star formation-the first genuine protogalaxies-is essential to 
understanding galaxy evolution. 

Pencil beam surveys have searched for z > 8 galaxies. In 
the deepest available near-infrared data — the Hubble Deep 
Fields North and South (HDF-N and - S), the Ultra Deep Field 
(UDF) and the UDF Parallel Fields— Bouwens et alJ (120051) 
find three possible z ~ 10 galaxies. In addition, observations 
of gravitationally lensed galaxies have also revealed a few po- 
tential z > 8 galaxies (Richard et al. 2006, 2008). However, 
these deep, small area surveys find faint candidates which 
are difficult to follow up, so none of the redshifts have been 
confirmed by measuring their spectral energy distributions 
(SEDs) at longer wavele ngths. This crucial te st has ruled out 
other z ~ 10 candidates (Bou wens et al.ll20p5h . including the 
one first presented in Dic kinson et all d2000l) . 

Our NICMOS Parallel Imaging Survey has the unique com- 
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bination of sensitivity and area to find rare, luminous J- 
dropout LBGs at 8 < z < 10. These candidates are typically 
#160 = 23 -25 (AB), so follow-up observations can confirm or 
rule out extremely high redshifts. In particular, these galax- 
ies can be observed with IRAC on Spitzer Space Telescope 
(hereafter Spitzer). We present here the first results of such a 
survey, where one possible z ~ 9 galaxy has been discovered. 
We use H = 70 kms" 1 Mpc" 1 , Q A = 0.7, and Q M = 0.3, and 
AB magnitudes are also used throughout. 

2. OBSERVATIONS & DATA REDUCTION 

We selected J-dropout candidates fro m the NICMOS Par- 
allel Imaging Survey ( ~ 135 arcmin 2 ; iTeolitz etalJfl 998: 
lYan et dJl2000t IColbertet al.ll200H iHenrv et alJl2007h . In 
this paper, we have adjusted the photometry by -0.16 and 
-0.04 mag in /no and H\w according to de Jong (2006). 
From more than seven thousand detected sources, four tar- 
gets were selected from the small subset (2%) with red /no - 
#160 (> 0.9) and no /no detection above 2a. 

We observed these NICMOS selected J-dropout candidates 
for four, five or six orbits each with NICMOS in the J\ 10 -band 
to confirm their r ed /no - i/160 color. The images were driz- 
zled according to Frucht er" & Hookl (119971) . and the resulting 
point-spread function has a FWHM^ 0.3". The 5 a sensitiv- 
ity, measured in a 0.8" diameter aperture is /no = 26.3- 26.5. 
The NICMOS colors were measured in apertures of this size. 

These four galaxies were also observed with the Infrared 
Array Camera (IRAC) on Spitzer, at 3.6, 4.5, 5.8, and 8.0 /im. 
The integration times were 0.9 ks pixel" 1 at 3.6 and 5.8/im, 
and 2.8 ks pixel" 1 at 4.5 and 8.0 /im. We used the Spitzer Sci- 
ence Center mosaicking software, MOPEX, with the drizzle 
algorithm to combine the images from the Basic Calibrated 
Data set. The output pixel size is 0.6", and resulting FWHM 
of the point-spread-function is 1.8" in the 3.6 and 4.5 /im 
bands, 2.0" in the 5.8 /im band, and 2.2" in the 8.0 /im band. 
The noise was measured by placing apertures in blank parts of 
the images. We chose a 2.4" diameter aperture to avoid con- 
tamination from nearby sources. For this small aperture, the 
corrections are 1 .8, 1 .9, 2.4, and 2.7 in the four IRAC bands 
(derived from point sources in our fields). The resulting 3 a 
sensitivities are 0.8, 0.8, 5.8, and 3.8 /iJy (total). 

3. RESULTS 
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Fig. 1 . — Postage stamp images of JD2325+1433 show that it is elongated in the NICMOS images. Images are 6" on a side. Note that the integration in 
/no is 10 times longer than in Hi^q, and the integration in IRAC1 (3.6 fim) is 13 times longer than in IRAC2 (4.5 fim), so the galaxy is indeed red in Jhq - Hi^q, 
and [3.6]-[4.5]. The coordinates of JD2325+1433 are 23h24m46.52s, 14d32m594s, bootstrapped to the 2MASS frame. 



The resulting photometry shows that three of the four 
sources have SEDs which are consistent with lower-z inter- 
lopers. One has /no - = 0.4, and it too blue to be a 
viable z > 8 candidate. Two other sources have rising SEDs 
into the IRAC bands, with H^o - [3.6] = 2-3. This makes 
them likely lower redshift interlopers, such as dusty starburst 
galaxies at z ~ 2-4 (e.g. IChen & Marzkell2004t [Labbe etal.1 
120051: iBrammer & van Dokkumll2007HStutz etalJl2008h . The 
most remarkable source, however, is JD2325+1433 which has 
Jno -#160 = 1-7, and a spectrum which is fiat (in f v ) through 
3.6 and 4.5 fim (much bluer th an the controversial z ~ 6.5 J- 
dropout candidate reported by Mo basher et al.ll2005l) . These 
are the colors that are expected for LBGs at 8 < z < 10. 

Initially, JD2325+1433 was only detected at 4.5 ^m (f v = 
1.29 ±0.29 fiiy), even though the 3.6 and 4.5 fim images 
reach the same sensitivity. Since this could be evidence for 
a Balmer break at a redshift of z ~ 9, we obtained deep (10 
hour) Spitzer Directors Discretionary Time (DDT) observa- 
tions at 3.6 and 5.8 fim to verify the break. The data are 
reduced in the same manner as described above, and the 3 
a sensitivity is 0.23 and 2.0 fiJy at 3.6 and 5.8 fim. The 
candidate, JD2325+1433 is detected at 3.6 fim, with f u = 
0.75 ±0.08 fiiy, confirming this break between 3.6 fim and 

4.5 fim. We do not detect JD2325+1433 at 5.8 fim, which is 
consistent with its flat spectrum. 

Postage stamp images of JD2325+1433 in Figure [T] show 
another galaxy approximately 2.5"to the east. To test whether 
there is significant contamination from this source, we per- 
formed a PSF subtraction to remove it. This assumes that 
the nearby source is unresolved the IRAC images, which 
is consistent with its NICMOS size (FWHM~ 0.5") and 
the observed IRAC FWHM (~ 1.8"). We then re-measure 
JD2325+1433 with the nearby source removed. The resulting 
flux density is f v = 0.82 and 1 .25 fiiy at 3.6 and 4.5 fim. This 
is consistent with the expectation that at most 2-3% of the 
light from a point source at the position of the nearby source 
would fall in the aperture that we used for JD2325+1433 — a 
fraction that corresponds to at most 0.05 fiiy. 

As an indep e ndent check, we used SExtractor 
dBertin & Arnoutsl [1996) to deblend the two sources. 
The resulting photometry, using the same 2.4 " diameter 
aperture, is consistent with aperture photometry when no 
correction for blending is made — /„ =0.81 and 1.33 fiiy at 

3.6 and 4.5 fim. Therefore, we conclude that contamination 
is negligible, and make no correction to the aperture photom- 
etry. We do, however, include the scatter derived from these 
three measurements in our IRAC error estimates, although 
it is only appreciable for the 3.6 fim photometry (0.04 fiiy, 
added in quadrature with 0.08 ^Jy). 



3.1. Matching IRAC & NICMOS Photometry 

To facilitate comparison with IRAC, we must estimate an 
aperture correction to account for missed flux in the NIC- 
MOS images. Although JD2325+1433 does not have the 
S/N required to characterize its morphology, we can place 
rough constraints on its size. We use GALFIT dPeng et all 
120021) to derive sizes and axial ratios for a few different as- 
sumed light profiles. For an exponential disk model, we find 
r s = 0.2 ± 0.03 " and a/b = 0.36 ± 0.09, and for a Gaussian 
we find a FWHM = 0.53 ± 0.06 ", and a/b = 0.36 ± 0. 1 1 . We 
also attempted to fit a de Vaucouleurs profile, but it was un- 
constrained. The uncertainties on these fits define the largest 
and smallest galaxy sizes, and therefore the largest and small- 
est aperture corrections that we can expect. To account for a 
possible de Vaucouleurs light profile, we add 0. 15 magnitudes 
of missed flux t o the largest expone ntial disk aperture correc- 
tion as found by Col bert et al.1 (120061) . We choose an aperture 
correction in the center of this range, and la uncertainties to 
span the entire range. For a 0.9" diameter aperture, this aper- 
ture correction is 1.18 ±0.12. We find that this aperture size 
minimizes the uncertainty on the total NICMOS flux density, 
which includes uncertainties for both the background and the 
aperture correction. We apply this aperture correction to Jno 
and //160, since the PSF is approximately the same for both. 

Finally, the total flux densities are: (Jno, H^o, [3.6], [4.5]) 
= (0.11, 0.50, 0.75, 1.29) ± ( 0.03, 0.09, 0.09, 0.29) fiiy. The 
3cr upper limits at 5.8 and 8.0 fim are 2.0 and 3.8 fiiy. 

3.2. On the Possibility of a Stellar Contaminant 

A stellar or substellar contaminant for JD2325+1433 is im- 
probable, for a several reasons. First, the source is likely 
resolved. In {sl3.ll we show that the semi-major axis has a 
FWHM ~ 0.5", while the NICMOS PSF has a FWHM of 
~ 0.35". Furthermore, the NICMOS images in Figure[TJshow 
that JD2325+1433 is elongated in the same direction in both 
the 7no and H\a) images. 

Even if JD2325+1433 were unresolved, its photometry can 
not be fit by a stellar contaminant. Of particular concern are 
brown dwarfs, which have molecular features that can imitate 
spectral breaks. While the reddest L-dwarfs can fit our NIC- 
MOS color, they can not simultaneously fit our IRAC pho- 
tometry (typically [3.6] - [4.5] < -0.5). Similarly, T-dwarfs 
can fit our IRAC photometry with [3.6] - [4.5] > 0.5 (AB), 
but with J - H < 0, they are much too blue to fit our NIC MOS 
observations dBurgasser et al.ll2006t iPatten et ai]|2006l) . This 
trend in the colors of late -type T-dwarfs is pr edicted to ex- 
tend to the coolest dwarfs (Baraffe et al. 2003; Barman, priv. 
comm.). Lastly, we note that other young stellar objects such 
as T-Tauri stars can not produce this red Jno - H\a) color (e.g. 
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FIG. 2. — Left- The photometry of JD2325+1433 shows two breaks in the spectrum, consistent with the Lyman and Balmer breaks and suggesting aredshift 
of z ~ 9. The black curve is the best fit model (z = 9.0 with an age of 64 Myrs, a constant SFR of 930 M Q yr 1 , A v = 1 .0, M s = -23.8, M = 5.4 X 10 10 M© 
and Z = O.2Z0). The grey curve is an alternate, poor fit at z = 2.9 (an instantaneous burst aged 360 Myrs, with Ay = 0.1, Mb = —21.1, M= 1.8 X 10 10 Mq and 
Z = Zq). Right- The results of our Monte Carlo simulation show probability as a function of redshift. More than 96% of solutions lie at z > 7. 

of JD2325+1433 is not well constrained, with values ranging 
from a few Myrs, to several hundred Myrs. The right panel 
shows a degeneracy between age and extinction. If we ex- 
clude models with ages less than 10 Myrs (since it is unlikely 
that all the stars in the galaxy can form within 10 Myrs), then 
we can also constrain the extinction to Ay < 1 . Lastly, we note 
that estimated stellar masses range from 10 10 to 10 11 M Q , for 
a Salpeter IMF. 

5. IMPLICATIONS 

If JD2325+1433 is at z ~ 9, its UV luminosity corresponds 
to 5L*(z = 3), and we can estimate a lower limit on the 
volume density of luminous 8 < z < 10 LBGs. 1 For this 
redshift range, the entire NICMOS Parallel Survey covers 
5.2 x 10 5 Mpc 3 (comoving) to Hk,q = 24. Beyond this, we 
are limited by /no - band sensitivity required to select sources 
with /no - Hi(,o > 0.9, as described in $2] At fainter magni- 
tudes (given the large inhomogeneity in pure parallel exposure 
times), 5 1 % of the total survey area has sufficient sensitivity 
to find galaxies like JD2325+1433. This implies that we have 
searched 2.6 x 10 5 Mpc 3 (comoving), and the resulting vol- 
ume density of 8 < z < 10 LBGs is > 3.8 x 10" 6 Mpc" 3 . 

The volume density implied here is rough l y cons istent with 
that at z ~ 3, as measured by ISteidel et all ( 119991) . This LF 
predicts (f> ~ 10" 6 Mpc -3 for galaxies with L > 5L*(z = 3), im- 
plying 0-0.75 such sources in our survey volume (assuming 
Poisson errors). On the other hand, z ~ 6 LF, as measured 
by Bouwens et al. (2006, 2007), predicts fewer luminous 
galaxies and therefore significantly smaller number for this 
survey (< 0.001). If this z ~ 6 LF is representative of z ~ 9 
LBGs, then JD2325+1433 is almost certainly an interloper at 
z ~ 2-3. However, this evolution in the LF is uncertain, and 
the z ~ 5 LF measured by Iwa ta et al.l (120071) shows no change 
in the number density of luminous galaxies from z ~ 3. Since 
the Iwata et al. survey contains four times more area than 
the Bouwens et al. survey, this measurement may be a more 
reliable probe of the most luminous sources. 

Indirect evidence also implies that z > 8 galaxies may be de- 
tectable in the NICMOS pure parallel imaging. At z ~ 5-6, 
the ages and masses of LBGs pr edict a more intense period 
of star formation at earli er times dYan et al.l l2006; Evl es et alj 
l2007t IStark et ai1l2007al) suggesting upward luminosity evo- 
lution at z > 6. This argument is also made for faint Lya 

1 Since we have not observed all NICMOS Parallel selected J-dropout 
galaxies with IRAC, several more galaxies like JD2325+1433 may exist in 
our survey. 



iHernandez et al. 2007), as their NIR colors are dominated by 
photospheric emission. 

4. CONSTRAINING THE REDSHIFT OF JD2325+1433 

We use the Hyp erz ph otometric redshif t code 
(iBolzonella et alJ 120001) with iBruzual & Charlotl (120031) 
stellar snythesis templates to estimate the redshift of 
JD2325+1433. We do not include the 5.8 /im limit, as no rea- 
sonable models predict a flux > 2 /xJy at 5.8 ^tm. Assuming 
a Salpeter initial mass function (IMF) and using the Padova 
1994 evolutionary tracks, we fit for red shift, varying age, 
mass normalization, extinction (Ay = 0-3; Calzett iet alj2000l 
law), and metallicity (Z = 0.005,0.02,0.2,0.4, 1.0 x Z Q ). 
We also allow three different star formation histories- an 
instantaneous burst, a constant star formation rate (SFR), and 
an exponentially declining SFR with an e-folding time of 100 
Myrs. With four data points, the problem is underconstrained, 
and we are not able to estimate metallicity or star formation 
history. However, we include these parameters so that we do 
not overlook some combination of parameters that is a good 
fit to our data at lower redshift (z ~ 2-4). 

In Figure [2] we compare the best fit model (black curve), 
which has z = 9.0, to the best possible fit that can be achieved 
for z < 6 (grey). This lower-z fit, which has z = 2.9, is sig- 
nificantly poorer, and can not reproduce the two breaks in the 
spectrum. To further test the significance of these fits, we gen- 
erated 10 5 Monte Carlo realizations of the photometry, simul- 
taneously varying the fluxes according to the measured un- 
certainties. We fit these data with Hyperz, as described above. 
The probability distribution in redshift space (see Figure [2] 
right), shows a peak at z = 9.0 and a much smaller peak at 
z ~ 2.5. More than 96% of solutions lie at z > 7. This is con- 
sistent with the observed 4.5 flux density lying only 1 J a 
away from the z = 2.9 model. For normally distributed uncer- 
tainties, the true 4.5 flux density will lie more than 1 .7 a 
away from our measurement 10% of the time. This implies 
that there is a 5% chance that the Balmer break is more than 
1.7 a larger than we have estimated, and 5% chance that there 
is no Balmer break at all. In the latter case, z ~ 2- 3 models 
will be a good fit to the photometry. 

In Figure [3] we show two dimensional confidence intervals 
in age/redshift space, and extinction/age. Estimated 1, 2, and 
3 a contours are determined from the probability at which 
68.3% (95.4%, 99.7%) of the Monte Carlo realizations have 
a higher likelihood. In the left panel we show that the z ~ 9 
solution is preferred at the 1 a level, while the 2 a interval in- 
cludes some z ~ 2.5 solutions. Both panels show that the age 
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FIG. 3. — Confidence intervals show that az ~ 9 solution is preferred at the Icr level, while age and extinction are not strongly constrained. Contours are 1, 2, 
and 3 a. On the left hand panel, contours are open because best-fit ages are restricted to be younger than the age of the universe. 



em itting galax i es. If any of the six candidates reported 
by IStark et al.1 ([2007b) are indeed at z > 8, then the space 
density of faint Lya emitters at this epoch must be larger 
than at z = 5.7 (as m easured by Malhotra & Rhoads l2004t 
Shimasaku et al. 2006). In conclusion, the density of galaxies 
and the evolution of the LF at z > 5 is sufficiently uncertain 
that the z ~ 9 interpretation of JD2325+1433 can not be ruled 
out. 

Constraining the numbers of galaxies at z > 8 is crucial to 
understanding the transition from a neutral to ionized inter- 
galactic medium (IGM). Currently, it is unclear if star forming 
galaxies at z ~ 6 are capable of reionizing the IGM, la rgely i n 
part due to uncertainties in the stellar IM F (e.g. ICharvl[2007l) . 
the escape fraction of ionizing photon s dMalkan et al.ll2003l : 
IShaplev et aLll2006h ISiana et alj|2007l) and the gas clumping 
factor. If JD2325+1433 is indeed at z ~ 9, the relatively large 
number of UV-luminous galaxies at this redshift will likely re- 



lax the requirement for strong evolution in these parameters. 

Finally, we note that further observations can resolve the 
degeneracy between the z ~ 9 and z ~ 2- 3 interpretations. 
The blue 7/ieo - IRAC colors constrain the SED so that we 
can predict Iab ^27-28 for the intermediate-z solutions. Of 
the 4% of Monte Carlo realizations (discussed in @ that have 
z < 7, none are fainter than I AB = 29, so deep optical imaging 
can definitively confirm the redshift of JD2325+1433. 
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